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PREFACE 

This  Technical  Memorandum  presents  an  overview  of  the  system  identification 
program  that  was  initiated  at  the  Naval  Air  Test  center  in  1970.  This  program  has 
resulted  in  an  advanced  system  identification  capability  that  has  a wide  spectrum 
of  applications  to  aircraft  flight  testing.  An  overview  of  maximum  likelihood 
parameter  identification  and  applications  to  aircraft  stability  and  control  flight 
testing  are  presented. 

This  Technical  Memorandum  was  prepared  for  presentation  to  the  AIAA 
Systems  and  Technology  Meeting  in  Dallas,  Texas,  on  27-29  September  1976.  The 
research  performed  to  develop  the  technology  reported  on  was  conducted  under  a 
series  of  AIRTASK  Assignments  and  Procurement  Requests  sponsored  by  the  Naval 
Air  Systems  Command  and  Office  of  Naval  Research,  respectively.  These 
programs  were  managed  by  Mr.  Ralph  A'Harrah  (AIR-53011)  at  the  Naval  Air 
Systems  Command  and  Mr.  Dave  Siegel  (ONR  Code  211)  at  the  Office  of  Naval 
Research.  A research  team  at  Systems  Control,  Incorporated,  headed  by  Dr.  W.  E. 
Hall,  developed  the  computer  program  SCIDNT. 
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Strike  Aircraft  Test  Directorate 
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The  development  of  system  identification  technology 
was  undertaken  to  provide  for  more  effective  aircraft 
flight  testing  oy  reducing  the  time  required  to  conduct 
specific  tests  and/cr  to  provide  for  a more 
comprehensive  data  analysis.  F-14A  and  TA-4J  flight 
test  results  presented  demonstrate  that  the  flight  time 
required  to  obtain  stability  and  control  data  can  be 
significantly  reduced  without  loss  in  accuracy  of 
conventional  flight  test  derived  parameters. 

Presentation  of  S-3A  and  EA-6B  system  identification 
results  demonstrate  that  this  technology  can  be 
successfully  used  to  update  the  aerodynamic  data  bases 
of  modern  jet  aircraft  from  flight  test  data.  These 
system  identification  results  are  compared  with  wind 
tunnel  data  and  flight  test  derived  parameters  to 
demonstrate  the  accuracy  of  this  new  technology. 
Applications  of  this  technology  to  integrate  several  areas 
of  aircraft  flight  testing  are  discussed. 
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L Primed  partial  derivative  of 

rolling  moment  with  respect 
to  x (where  x is  p,  r,  , 

£ , Sa,  8 K , 8 sp)  sec 

M Partial  derivative  of  pitching 

moment  with  respect  to  x 
where  x is  q,  a , u , 8e  ) sec 

N Primed  partial  derivative  of 

yawing  moment  with  respect 
to  x (where  x is  p,  r,/3  , , 

8<j,Sr,  8Sp)  sec 

X Partial  derivative  of  X force 

with  respect  to  x (where  x is  _j 

q,  a , u ,®«  ) sec 

Y Partial  derivative  of  sideforce 

with  respect  to  x (where  x is  _j 

p,  r,/3  ,80.8  R,8Sp)  sec 

Z Partial  derivative  of  Z force 

with  respect  to  x (where  x is 

q,  <*  » u .8  e ) 

C^  Non-dimensional  partial  deriv- 

‘x  ative  of  K force  or  moment  (L, 

M,  N,  X,  Y,  Z)  with  respect 
to  state  vector  x or  control 
vector  8 ) 

Subscripts 

m Measured 

T True 

( L Trim  condition 

o 

Superscripts 

T Transpose 

Background 

The  Naval  Air  Test  Center  (NATO  initiated  a 
program  in  1971  to  develop  advanced  system 

identification  techniques  for  use  in  flight  testing 
aircraft.  This  program  has  been  a coordinated  effort 
among  the  Naval  Air  Test  Center,  Naval  Air  Systems 
Command,  Office  of  Naval  Research,  and  Systems 
Control,  Incorporated.  The  major  objective  of  the 
development  of  this  technology  has  been  to  provide  for 
mor<1  effective  flight  testing  by: 

a.  Improving  safety  of  flight. 

b.  Reducing  cost  and/or  time  associated  with  design, 
flight  tests,  and  certification  of  aircraft. 

c.  Improving  data  analysis  accuracy. 

d.  Providing  a basis  for  more  comprehensive  analysis 
of  mission  suitability. 

e.  Providing  for  the  acquisition  of  flight  test  data 
for  use  in  flight  trainers. 

f.  Providing  an  accurate  data  base  for  the  design 
and/or  modification  of  advanced  flight  control  systems. 

To  date,  the  application  of  system  identification  to 
Navy  flight  testing  has  been  primarily  in  the  areas  of 
stability  and  control  testing  and  the  updating  of  aircraft 
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aerodynamic  data  bases  for  use  in  system  redesign  or 
modification.  For  example,  the  determination  of  the 
compliance  of  an  aircraft  stability  and  control 
characteristic  with  the  requirements  of  Military 
Specification  MIL-F-8785B  is  a costly  and 

time-consuming  facet  of  aircraft  flight  testing. 
However,  a considerable  portion  of  the  stability  and 
control  flight  program  can  be  eliminated  through  use  of 
system  identification  by  using  this  technology  to  extract 
the  aerodynamic  stability  and  control  derivatives  from  a 
limited  number  of  flight  tests.  These  stability  and 

control  derivatives  are  then  used  to  verify  the  aircraft's 
compliance  with  the  Military  Specification  requirements. 

Overview  of  NATC  System  Identification  Program 

In  order  to  accomplish  the  overall  objectives  of  the 
development  of  this  technology  (items  a.  through  c.  of  the 
preceding  section),  a comprehensive  system  identification 
capability  will  be  required.  In  general,  a total  system 
identification  capability  is  considered  to  consist  of: 

a.  Design  of  experiments  (input  design). 

b.  Model  structure  determination. 

c.  Parameter  identification. 

The  algorithms  that  NATC  is  developing  to  formulate 
this  system  identification  capability  are: 

a.  Linear  maximum  likelihood  parameter  identifica- 
tion. 

b.  Nonlinear  maximum  likelihood  parameter  identifica- 
tion. 

c.  Data  consistency  (state  estimation). 

d.  Instrumental  variable  parameter  identification 
(real-time  parameter  identification). 

e.  Transfer  function  analysis. 

f.  Time  series  analysis. 

g.  Vehicle  dynamics  simulation. 

h.  Model  structure  determination. 

The  one  key  element  in  the  formulation  of  all  of  these 
algorithms  is  that  they  have  been  programmed  in  a 
general  format  such  that  they  can  be  easily  modified  for 
the  analysis  of  any  type  of  system.  Item  a.  is  an 
operational  program,  and  items  b.  through  f.  have  been 
installed  on  the  NATC  computer  system  and  are  currently 
under  evaluation.  Items  g.  and  h.  are  currently  under 
development  and  are  planned  for  completion  by  the  end  of 
this  year. 

Parameter  Identification  Procedure  Used  in  this  Analysis 

The  results  presented  in  this  paper  are  based  on  a 
linear  maximum  likelihood  parameter  identification 
computer  algorithm  (SCIDNT  HD  currently  being  used  in 
the  analysis  of  flight  test  data  at  the  Naval  Air  Test 

Center.  This  program  provides  for  the  estimation  of: 

a.  Coefficients  of  linear  systems. 

b.  Instrumentation  system  scale  factor  errors,  biases, 
and  lags. 

c.  Gust  time  history  characteristics. 
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Thus,  if  a system  is  modeled  in  conventional  state-space 
notation: 

State  Equation 

« = F»  ♦ G8 ♦I'**  (1) 

where: 

x is  a nxl  state  vector 

8 is  a I xl  control  vector 

w is  a qxl  random  process  noise  vector 

F is  a nxn  matrix  of  stability  derivatives 

G is  a nx  / matrix  of  control  derivatives  and 
corrections  for  initial  out-of-trim  conditions 

r is  a nxq  process  noise  distribution  matrix 

Measurement  Equation 

y = Hx«-D$*b*v  (2) 

where: 

y is  a rxl  measurement  vector 

B is  a rxl  measurement  bias  vector 

v is  a rxl  measurement  random  noise  error  vector 

H is  a rxn  matrix  relating  the  measurements  to  the 
state  vector 

D is  a rx  l matrix  relating  the  measurements  to  the 
control  vector 


and 


n is  the  number  of  states 

t is  the  number  of  controls 

r is  the  number  of  measurements 

and  w and  v are  Gaussian  random  noises  which  have  zero 
mean,  are  uncorrelated,  and  have  power  spectral 
densities  Q and  R,  respectively. 

Then,  this  parameter  identification  algorithm 
provides  for  estimates  of  the  elements  which  make  up  F, 
G,  H,  D,  and  B.  In  addition,  estimates  are  made  of  the 
power  spectral  density  of  the  process  noise  w and  the 
measurement  noise  v.  This  identification  algorithm  is 
programmed  in  a general  format  such  that  the  state 
equation  and  measurement  equation  can  be  easily 
modified.  Thus,  this  program  is  easily  changed  to 
account  for  any  differences  in  aerodynamics  or  control 
system  characteristics  between  airplanes  or  to  model  any 
type  of  system  such  as  an  aircraft's  power  plant. 

In  general,  system  identification  is  the  process  of 
estimating  from  a given  set  of  input/output  data  either 
the  required  structure  of  a system  model  or  the 
parameters  of  a prespecified  model.  The  general  process 
of  parameter  identification  is  illustrated  in  figures  1 and 
2.  Figure  1 is  an  overview  of  the  elements  in  parameter 
identification,  and  figure  2 is  a schematic  of  the 
iterative  nature  of  the  algorithm.  As  shown,  the 
parameter  estimation  procedure  begins  with  the  system 
or  aircraft  being  disturbed  from  some  initial  condition  by 
a pilot  input.  The  aircraft  response  to  this  input  is  then 


recorded  on  the  instrumentation  system  (in  general,  this 
measurement  can  be  corrupted  by  external  disturbances 
and/or  instrumentation  noise).  This  measured  flight 
response  is  compared  with  a computed  response  based  on 
a mathematical  model  of  the  aircraft  to  form  a response 
error.  This  response  error  is  formulated  into  a criterion 
function  which  is  used  in  an  estimation  algorithm  to 
update  the  initial  guess  of  the  parameter  values.  This 
procedure  is  repeated  until  the  response  error  ideally 
goes  to  zero,  at  which  point,  the  'best'  estimate  of  the 
aircraft  parameters  is  obtained. 

In  the  maximum  likelihood  algorithm  used,  the 
criterion  function  is  the  likelihood  function 


(3) 


which  is  defined  as  the  conditional  probability  of  the 
measurements  (yi  having  occurred  given  the  parameter 
set  9 ( Q is  the  set  of  parameters  defined  from  the 
system  state  and  included  in  the  feasible  set  0 ).  Thus, 
for  any  given  set  of  values  of  the  parameter  6 , we  can 

assign  a probability  of  p [^/q J to  each  outcome  y.  If 

the  outcome  of  an  actual  flight  test  is  y,  it  is  of  interest 
to  know  which  set  of  values  might  have  led  to  these 
observations.  In  order  to  accomplish  this,  the  maximum 
likelihood  method  finds  a s**t  of  parameters  to  maximize 
the  likelihood  function  such  that 


(4) 


0 s max  P \/9) 

0«0 

In  other  words,  the^ probability  of  the  outcome  y is  higher 
with  parameters  9 in  the  model  than  with  any  other 
values  of  parameters  from  the  feasible  set.  The 
likelihood  function  used  is  given  by: 


‘ (y(t|)/y(t,-i),e)  = 


exp{-^-«T(i)BH(.)w(i)} 
(2xr)M/2 1 B(.)  I'/2 


or  in  more  convenient  log  form: 


where 


°(i)=y(i)-y(i) 


(5) 


(6) 


(7) 


and  m is  the  number  of  measurements  and  B 
approximated  by 


N 


8»rr  £ w(i)w  (i) 


N 


i=l 


(8) 


The  optimization  procedure  used  to  find  the  parameter 
set  is  the  modified  Newton  Raphson  technique 


where 


A A 
Sj.i'Sj.Ae 


A0=-  | 


d2log  P(y/e)|~'dlog  P(y /8) 

M ( de, 


(9) 


(10) 


and  i is  the  ith  iteration  and  j and  k denote  the  jth  and 


kth  parameters. 


(3) 
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Applied  Huns  of  System  Ulentifn  ation  to  Determining 
Stability  and  Control  Flight  Test  Requirements 

The  military  flying  qualities  specification, 

MIL-F-8785B,  formulates  many  stability  and  control 
requirements  in  terms  of  transfer  function 

characteristics  and  mode  ratios.  For  example,  many 
dynamic  longitudinal  and  lateral  directional 
characteristics  are  specified  in  terms  of  transfer 
function  numerator  (zero)  and  denominator  (pole)  root 
locations  (frequency,  damping,  and  time  constants). 
Static  stability  characteristics  such  as  longitudinal 
control  (elevator  position)  variation  with  airspeed  are 
related  to  the  appropriate  transfer  function  evaluated  at 
steady  state  conditions  (the  Laplace  operator  s = 0). 
Specification  requirements  for  the  ratios  of  normal 
acceleration  to  angle  of  attack  (n*/a  ) and  roll  angle  to 

sideslip  angle  ( yfei  can  be  determined  using  mode  ratios 
evaluated  at  the  correct  frequency.  Compliance  with 
these  types  of  specification  requirements  can  easily  be 
determined  through  use  jf  system  identification* 

If  a simplified  form  of  equations  (1)  and  (2)  is 
considered 


and  the  denominator  polynominal  (characteristic  equation) 
is 

A=(s2*2;pwnps^np2)(s2’2Up“ns^wnsp^)  (lg) 

Phugoid  Mode  Short  Period  Mode 


Similarly,  the  normal  acceleration  to  elevator  input 
transfer  function  is  given  by: 


nZ 

nz(s)  NSe 

I^sTT 


(19) 


where  the  numerator  polynominal  is  of  the  form 

Ns2=KQj,s(s.Vr0zi)(s.l/raz2)(s.l/rQ,J 


(20) 


Thus,  the  dynamic  longitudinal  stability  specification 
requirements  that  can  be  determined  by  evaluating 
equations  (16)  through  (20)  are  the  phugoid  frequency 

(WnP  ) afld  damping  (Cp)  and  short  period  frequency 
(u'nsp  ) and  damping  (C$p)/  This  is  further  illustrated  by 


x = Fx ♦ G8  (11) 

y = Hx^D8  (12) 

The  Laplace  transform  solution  of  equation  (11)  for  zero 
initial  conditions  is  given  by: 

x(s)*(sI-FflG8(s)  (13) 

where  I is  the  identity  matrix.  This  expression  is  now 

substituted  into  the  Laplace  transform  of  equation  (12) 
to  obtain 

y(s)*H(»I-F)'G8(s)*D8(s)  <14> 

A general  expression  for  the  output  to  input  transfer 

functions  of  the  system  described  by  equations  (11)  and 
(12)  is  then  given  by: 

-^=H(sI-F)_lG.D  (15) 


If  the  elements  in  the  matrices  F,  G,  H,  and  D can  be 
estimated  from  flight  test  data  using  system 

identification  techniques,  it  is  easily  seen  that  equation 
(15)  can  be  used  to  determine  the  stability  and  control 
specification  requirements  previously  mentioned.  (It 
should  be  noted  that  this  equation  is  solved  on  a digital 
computer  to  obtain  numerator  and  denominator  roots  and 
classical  approximations  for  these  terms  are  not  made.) 
Specific  examples  for  determining  longitudinal  and 
lateral-directional  dynamic  and  static  stability 
characteristics  are  now  presented. 

Longitudinal  Specification  Requirements 

The  angle  of  attack  to  elevator  input  transfer 
function  is  given  by: 

•(«)  % 

SSrt'T  (16> 


where  the  numerator  polynominal  is  of  the  form 

NV«a<s.  lV(s2.2{a“na*.«na) 


(17) 


considering  a classical  approximation  to  these  terms,  from 


which  it  is  seen  that  once  the  stability 
derivatives  are  known,  it  is  a simple  matter 
these  specification  requirements. 

and  control 
to  calculate 

r - 1 l^tXa-g)  , 

tp  2-np)  Ma 

(21) 

V«9<MUaZa-Zu) 

and 

(22) 

2<*  <Z„.Mq) 

nsp 

(23) 

"nSp2«Z0Mq-Ma(Zq.|) 

(24) 

The  normal  acceleration  sensitivity  specification 
requirement  can  be  obtained  by  forming  the  mode  ratio 

( Z/a  ) from  equations  (16)  and  (19)  and  evaluating  it  at 
the  short  period  root  location: 


Static  longitudinal  stability  requirements  in  terms  of 
elevator  position  gradients  with  airspeed  can  be 
determined  by  evaluating  the  airspeed  ( u ) to  elevator 
input  transfer  function  at  steady  state  conditions. 

u(s)  _ ^8e  . l4u(S*l/Tu|)(S*l/TUg)  (26) 

8e(*f  fl  A 


Evaluating  (26)  at  s = 0 

u(s) 

».<•> 

s*o 

which  is  equivalent  to 


KuC/TU.H'/TU;) 

“"•p2“"o* 


u(s)  -(-ZoM8e+MaZ8e) 

Vs>  s=p  ” <MaZu-MuZa) 


(27) 


(28) 
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Thus,  equation  (28)  can  be  interpreted  as  the  gradient  of  Flight  Test  Results 

airspeed  with  elevator  position  during  a static 

longitudinal  test.  Reductions  in  Flight  Testing  Using  System  Identification 


Lateral-Directional  Specification  Requirements 


In  the  lateral-directional  axes,  MIL-F-8785B  sets 
new  requirements  for  Dutch  roll  damping  (‘"fid  I and 

frequency  (C<j),  spiral  mode  (l/rs),  roll  mode  (1/rq), 

and  roll  rate.  In  addition  to  the  updated  military 
specification  requirements  in  the  lateral-directional 
axes,  there  are  also  new  parametric  requirements  in  th** 
detail  specification  for  the  S-3A  and  F-14A  airplanes. 
These  new  requirements  are  in  the  form  of  the  Dutch  roll 


coupling  parameter 
parameter  j 


and  the  Dutch  roll  excitation 


These  new  specification 


requirements  in  the  lateral-directional  axes  are  difficult 
to  determine  accurately  because  the  effects  of  the 
spiral,  roll,  and  Dutch  roll  modes  cannot  be  easily 
separated  using  ccnventional  data  techniques.  However, 
these  new  requirements  can  easily  be  determined  if  the 
roll  rate  to  aileron  input  transfer  function  is  evaluated 
using  estimated  stability  and  control  derivatives. 


,<s> 


8a*s)  (s.'/TpMs.'/rjHs^SC^nd5*"^ 


(29) 


Compliance  with  the  specification  requirements  can  be 
determined  from  the  estimated  transfer  function 

parameters  “nd  , , VT  , and  - . This  is 

R *s, 

illustrated  by  considering  a classical  approximation  to 
the  roll  mode  time  constant 


(30) 


In  order  to  determine 


«d 

Kss 


from  these  data,  the  matched 


transfer  function  poles  and  zeros  are  plotted  on  a s-plane 
as  shown  in  figure  3.  The  terra  is  then  determined  as 

the  residue  measured  from  the  Dutch  roll  pole  and  is 
given  by: 


K 


d 


Dutch  = 

Roll 

Residue 


0 b 

ewnd“d 


(31) 


where  a,  b,  and  e are  defined  in  figure  3.  The  term 

is  the  steady  state  residue  and  is  measured  from  the 
origin,  assuming  that  l/r#  = 0. 


K 


ss 


Steady 

State 

Residue 


(32) 


Thus,  K^/Kw  is  determined  as  the  ratio  of  equations  (31) 
and  (32)* 


System  identification  technology  can  be  used  to  reduce 
the  flight  time  required  to  obtain  data  for  determining 
compliance  with  stability  and  control  specification 
requirements.  This  is  accomplished  by  using  this 
technology  to  extract  from  a limited  number  of  flight 
tests  the  stability  and  control  derivatives  which  are  then 
used  to  determine  compliance  with  specification 
requirements  which  normally  require  multiple  tests  at 
each  trim  flight  condition.  In  order  to  accomplish  this 
objective,  it  was  necessary  to  conduct  a flight  program  to 
determine  the  optimal  flight  inputs  for  parameter 
(4) 

identification  analysis.  The  result  of  this  research 
effort  and  additional  follow-on  work  has  been  the 
selection  of  a sequential  aileron-rudder  doublet  for 
lateral-directional  stability  analysis  and  an  elevator 
doublet  or  sine  wave  for  longitudinal  short  period  analysis 
(a  sine  wave  appears  to  have  no  advantages  over  a doublet 
input  for  linear  analysis).  Parameter  identification 
analysis  of  the  phugoid  mode  requires  use  of  the 
conventional  test  technique. 

The  use  of  these  'optimal'  inputs  to  reduce  flight  tests  is 
demonstrated  in  table  I.  As  shown,  system  identification 
techniques  can  be  used  to  reduce  the  total  number  of 
maneuvers  required  in  determining  compliance  with  the 
longitudinal  arid  lateral-directional  stability  and  control 
specification  requirements  previously  noted  by  a factor  of 
3 at  each  individual  trim  point  (a  reduction  from  9 
maneuvers  to  3 maneuvers).  This  reduction  in  the  number 
of  maneuvers  required  results  in  a savings  of 
approximately  75  percent  in  flight  test  time  (based  on 
TA-4J  flight  tests  conducted  by  the  U.S.  Naval  Test  Pilot 
School).  The  TA-4J  was  then  used  in  a test  designed  to 
provide  flight  data  for  determining  the  aircraft's 
characteristics  over  an  airspeed  range  at  one  altitude. 
This  consisted  of  accelerating  the  aircraft  from  165  knots 
indicated  airspeed  to  the  maximum  level  flight  airspeed 
(at  15,000  feet)  and  collecting  data  at  50  knot  increments. 
Tests  were  conducted  at  seven  specific  trim  points  as 
shown  in  table  2.  At  each  test  point,  the  aircraft  was 
stabilized  (trimmed)  and  an  aileron -rudder  sequential 
doublet  and  elevator  doublet  inputs  were  made  by  the 
pilot.  (These  tests  did  not  include  phugoid  maneuvers.) 

The  total  flight  test  time  required  for  these  tests  was 
13  minutes.  Using  conventional  flight  test  procedures,  the 
time  required  to  obtain  the  flight  data  for  determining  the 
same  specification  requirements  is  estimated  to  be  112 
minutes.  Thus,  if  phugoid  test  data  are  not  required,  the 
reductions  in  flight  time  are  even  more  dramatic.  This 
saving  in  flight  time  is  due  to  the  reduction  in  the  time 
required  to  conduct  the  maneuvers  and  the  time  required 
to  establish  a fewer  number  of  precise  trim  points. 

A similar  set  of  tests  was  conducted  during  the 

Technical  Evaluation  of  the  F-14A  airplane. ^ These 
tests  are  summarized  in  table  3.  As  previously 
demonstrated,  a significant  reduction  in  flight  test  time 
was  achieved  using  the  system  identification  approach  to 
obtaining  flight  data  (a  reduction  of  32  to  8 maneuvers 
was  achieved  during  these  tests).  Total  test  time  was  8 
minutes. 

To  demonstrate  the  accuracy  of  this  approach  to 
flight  testing,  these  F-14A  data  were  used  to  estimate 
the  longitudinal  short  period  characteristics.  The  body 
fixed  axis  system  equations  used  in  this  analysis  are 
presented  below.  Elements  of  the  state  equation  are: 
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Slats  Vector  „=  [o,u,q,fi]T 
Control  Vector  8 = [s e , |]  T 
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(34) 
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The  measurement  equations  are*. 

Measurement  Vector  y=  [am,um,qm  ,Sm,0Zm]T  (38) 

(39) 


Bias  Vector 


b = [ba.bu'bq.be.baz]' 


Random  Measurement 
Error  Vector 


Vz  j“a'uu'l,q  >v8  ,uaz]T 


(40) 


Figure  1 presents  the  F-14A  stability  and  control 
derivative  estimates  and  associated  estimate  confidence 
bound  . Figure  5 presents  the  short  period  characteristics 
computed  from  these  derivative  estimates  and  compares 
them  with  conventional  flight  test  results.  As  shown, 
excellent  agreement  is  obtained  between  the  conventional 
results  and  the  parameter  identification  results 
(conventional  results  are  based  on  classical  hand 
measurement  techniques).  There  is  a significant  point  to 
make  here  in  that  there  is  no  comparison  of  conventional 
flight  test  determined  n*/a  with  the  parameter 
identification  result  because  it  would  have  required 
another  test  to  obtain  the  conventional  data.  A sample 
time  history  match  is  presented  in  figure  6* 

The  examples  presented  demonstrate  the  accuracy  and 
feasibility  of  using  system  identification  techniques  to 
improve  the  efficiency  of  stability  and  control  flight 
testing.  Of  course,  the  examples  presented  here  represent 
only  a portion  of  the  stability  and  control  testing 
requirements  and  thus  savings  in  flight  time  will  not  be  as 
dramatic  when  considering  the  total  testing  requirements. 
A previous  survey  of  Navy  stability  and  control  flight  test 
programs  conducted  at  NATO  indicated  that  70  - 90 
percent  of  the  testing  in  large  scale  programs  is  devoted 

to  specification  testing.  Using  system  identification,  it 
is  estimated  that  20-30  percent  of  this  portion  of  the 
flight  program  could  be  eliminated. 
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Verification  of  S-3A  Power  Approach  Characteristics 


(41) 


(42) 

(43) 

(44) 


In  order  to  improve  the  carrier  suitability  of  the  S-3A 
airplane,  a program  was  initiated  by  the  Naval  Air 
Systems  Command  to  identify  the  origins  of  any  S-3A 
carrier  approach  difficulties  and  to  solve  them.  A portion 
of  this  program  was  to  apply  the  advanced  system 
identification  techniques  being  developed  at  NATC  to 
S-3A  power  approach  flight  test  data  and  compare  the 
resulting  characteristics  with  aerodynamic  data  the 
airframe  contractor  is  using  to  describe  the  airplane.  This 
verification  of  the  S-3A  data  base  was  accomplished  by 
comparing  original  contractor  data  with  NATC  and 
airframe  contractor  parameter  identification  results. 
(The  airframe  contractor  was  conducting  a parameter 
identification  analysis  concurrent  with  the  independent 

NATC  effort-^)  Results  from  the  NATC  analysis  have 

(8) 

been  published  and  are  summarized  below. 


(45)  state  space  model  used  in  this  analysis  takes  on  the 

following  form: 


(46) 

(47) 


Equation  of  Motion: 

State  Vector  x = [p,rj3,</>]T 
Control  Vector  8 = [8a,8*.l,8sp]T 


(49) 

(50) 


A general  form  of  the  accelerometer  measurement 
equation  (fifth  row  of  the  H and  D matrices)  is  given  by: 

°zm  = °zT',xq  + ba2*''az  <«) 

Thus,  using  these  equations  in  the  identification 
algorithm,  this  set  of  data  was  analyzed  assuming  that 
the  short  period  and  phugoid  modes  are  uncoupled 
(X  u = X a = Xq=  X g e = Z u = M u = 0).  The  results 

from  these  tests  are  presented  in  figures  4 and  5. 
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where  L , N , Y , and  are  initial  conditions  and 

o o o o 

primed  derivatives  are  defined  as: 
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and: 
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where  the  subscript 

(x)  denotes  p,  r, /3,S0  , 

and 

D(5.3)  = u0Y0-lzL'oe,No 

(67) 

^sp  . In  addition,  an  aerodynamic  spoiler  lag 

applied  to  spoiler  measurements. 

was 

0(5,4)  = u.Y*  -1,  L‘s  il.Ng 
’ 0 ®sp  2 °sp  * ®Sp 

(68) 

*sPlag=  ■ 

©Vag^pM 

(55) 

A general  form  of  the  lateral 
measurement  equation  can  be  written  as: 

acceleration 

where: 

0ym=0yT"p|2  *rl*  +b°/u°y 

(69) 

vTrue  (m  knots) 

This  correction  is  required,  especially  at  low  speeds, 
because  of  the  time  required  for  lift  to  build  up  after  a 
spoiler  input  (aerodynamic  lag).  The  solution  of  equation 
(55)  ( ®*P(_AG  ^ was  use(*  as  the  input  for  spoiler 

(equation  (50))- 
Measurement  Equations: 


Parameter  identification  analysis  results  are  presented 
in  table  4 and  figures  7 and  8.  As  shown,  final  airframe 
contractor  and  NATC  identification  results  are  in  general 
agreement  but  show  some  significant  differences  with  the 
original  data.  Based  on  this  comparison  and  the  time 
history  matches  presented  in  figure  8,  it  was  concluded 
that  the  original  data  base  did  not  accurately  represent 
the  lateral-directional  aerodynamics  of  the  S-3A.  Thus, 
the  NATC  parameter  identification  analysis  essentially 
verifies  the  final  airframe  contractor  aerodynamic  data 
base  of  the  S-3A  (excluding  spoiler  sideforce 
characteristics). 


Measurement  Vector  y = [Pm.  rm  ,0m  ,0^] T ^ 

g^-^t0r  b = [bp,br,b0,VbOy]T  <58) 
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As  with  the  use  of  any  new  technology,  additional 
comparisons  are  desirable  for  obtaining  confidence  in 
experimental  results.  These  comparisons  are  presented  in 
table  5 and  figures  9 and  10.  Table  5 presents  a 
comparison  of  lateral-directional  modal  characteristics 
(stability  and  control  specification  requirements).  Figure 
9 presents  a comparison  of  parameter  identification 
results  and  parameter  estimates  based  on  a flight  test 
differential  engine  thrust  sideslip  maneuver.  Figure  10  is 
a comparison  of  flight  test  results  for  a steady  heading 
sideslip  maneuver  (lateral-directional  static  stability 
characteristics)  and  parameter  identification  results. 
These  three  comparisons  show  good  agreement  between 
the  parameter  identification  and  conventional  flight  test 
results  and  in  this  light  verify  the  parameter 

identification  estimates. 


(h(5,I)  H(5,2)  H(5,3) 

oj 

EA-6B  Catapult  Launch  Capabilities 

where: 

H (5,1)  = u0Yp  -lz  L'p  ♦!„ 

N' 

(61) 

A program  was  initiated  by  the  Naval  Air  Systems 
Command  to  determine  EA-6B  catapult  launch 

r 

capabilities  as  limited  by  an  engine  failure  immediately 

H(5,2)  = u0Yr-lz  L‘r  +lx 

N’r 

(62) 

following  launch.  Part  of  this  program  was  to  determine 
from  flight  tests  the  longitudinal  aerodynamic 
characteristics  of  this  airplane  for  use  in  a catapult 
launch  simulation.  This  simulation  will  be  used  to  define 
critical  flight  areas  prior  to  any  actual  flight  tests.  The 

H(5,3)  = u0Y^g -lz  L^g  + I, 

(63) 

results  from  this  program  are  to  be  published  and  are 

(o) 

summarized  below. 
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The  slat**  equations  used  in  the  analysis  of  the  F-14 
airplane  were  augmented  in  this  case  to  include  a lag  on 
the  angle  of  attack  vane. 

v-ejs-ftMtxii). 


The  measurement  equation  for  the  angle  of  attack 
vane  was  then  modified  as  f- Hows: 


Jm=Ka<VtVna 


Parameter  identification  estimates  for  the  stability  and 
control  derivatives  and  comparisons  with  conventional 
flight  test  results  and  wind  tunnel  data  are  presented  in 
figure  11  and  table  6.  As  shown  in  figure  11,  reasonable 
agreement  between  the  Z force  derivatives  was  obtained; 
however,  there  is  a large  difference  in  the  estimates  for 
the  pitching  moment  derivatives.  The  largest  difference 
vs  vn  the  pitch  damping  derivative  ( Cmq^Cmg ) which 

would  indicate  a significantly  higher  level  of  short  period 
damping  than  is  obtained  from  the  wind  tunnel  data. 
(This  is  specifically  pointed  out  in  the  following 
paragriphs.)  For  the  phugoid  mode,  close  agreement  is 
obtained  between  the  parameter  identification  estimate 
and  conventional  flight  test  results  for  the  parameters 
C*u  and  CZy  . Representative  time  history 

comparisons  for  short  period  and  phugoid  aircraft 
responses  are  presented  in  figures  12  and  13. 

Data  are  now  presented  to  demonstrate  that  these 
stability  and  control  derivative  estimates  can  be  used  to 
determine  compliance  with  dynamic  and  static 
longitudinal  stability  specification  requirements. 

Comparisons  of  dynamic  stability  characteristics  are 
presented  in  table  7 and  show  good  agreement  between 
parameter  identification  and  conventional  flight  test 
results  for  both  the  phugoid  and  short  period  modes. 
However,  comparison  with  wind  tunnel  data  shows  that 
short  period  damping  estimates  are  approximately 
one-half  of  the  parameter  identification  and 

conventional  flight  test  results*  As  previously 

discussed,  this  would  indicate  that  the  wind  tunnel  value 
for  the  pitch  damping  derivative  is  low.  Comparisons  of 
static  longitudinal  characteristics  in  figure  14  show  good 
agreement  between  conventional  results  and  the  elevator 
to  airspeed  gradient  computed  using  the  parameter 

identification  stability  derivative  estimates*  ^ 
Integrated  Flight  Testing 

This  paper  has  dealt  with  the  savings  that  can  be 
realized  in  stability  and  control  flight  testing  by  utilizing 
system  identification  technology;  however,  an  even  more 
dramatic  reduction  In  flight  test  costs  and  time  could  be 
achieved  by  an  integrated  approach  to  flight  testing. 
This  integrated  flight  testing  will  be  made  possible  by  a 
further  growth  in  data  analysis  technologies,  such  as 
system  identification  and  dynamic  performance.  This 
can  be  readily  illustrated  by  considering  Navy  vehicle 
dynamics  tests  which  consist  of: 

a.  Aerodynamics. 

b.  Stability  and  control. 

c.  Performance. 

d.  Automatic  flight  control  system. 


e.  Automatic  carrier  landing  system* 

f.  Propulsion. 

g.  Structures  and  flutter. 

Under  current  test  philosophy,  these  tests  require  at  least 
four  airplanes  devoted  to  the  flight  tests  development 
program.  Considering  the  state  of  current  data  analysis 
technology  in  system  identification  and  dynamic 
performance,  it  is  possible  to  integrate  these  test 
requirements  to  form  a reduced  flight  program  (i.e.,  a 
reduction  in  both  flight  tests  and  required  aircraft  could 
be  achieved).  This  is  illustrated  by  considering  the  flight 
profile  in  figure  i5  which  consists  of  primarily  a series  of 
level  flight  accelerations/decelerations  and  constant  Mach 

climbs  and  descents-  The  aircraft  is  stabilized  at 

various  points  during  the  acceleration/deceleration 
maneuvers  and  various  maneuvers  are  performed  to 
collect  data  for  military  specification  compliance 
determination.  For  example,  if  these  maneuvers  consisted 
of  an  elevator  doublet,  aileron-rudder  double*  phugoid, 
and  engine  Bodie,  then  the  majority  of  stability  and 
control,  performance,  propulsion,  and  automatic 
stabilization  requirements  could  be  determined  using 
advanced  data  analysis  techniques.  (The 

acceleration/deceleration  data  are  used  in  dynamic 
performance  methods  to  estimate  aircraft  performance 
characteristics.)  Although  the  examples  used  do  not  cover 
all  of  the  test  requirements  for  these  areas,  such  as 
specialized  automatic  control  functions,  the  tests 
conducted  in  this  flight  profile  could  be  easily  modified  to 
include  them.  (This  would  extend  the  number  of  flights 
required  to  complete  the  flight  profile  for  one  aircraft 
loading  and  configuration.) 

Concluding  Remarks 

The  use  of  system  identification  technology  to  provide 
for  more  effective  aircraft  stability  and  control  flight 
testing  has  been  demonstrated.  This  is  accomplished  by 
either  improving  the  efficiency  of  the  flight  test  and/or 
providing  for  a more  comprehensive  data  analysis.  Thus, 
the  application  of  this  technology  to  flight  testing 
provides  for  an  in-depth  understanding  of  the  cause  and 
effect  relationship  in  aircraft  stability  and  control 
characteristics.  The  one  remaining  objective  to  be 
reached  in  this  program  is  to  apply  this  technology  in  a 
large  scale  flight  program  to  update  the  aerodynamic  data 
base  of  am  aircraft  throughout  its  flight  envelope.  To 
meet  this  objective,  current  plans  call  for  the  application 
of  this  technology  in  the  development  programs  for 
automatic  carrier  landing  systems,  TA-4KU  aircraft, 
HARRIER,  and  the  Navy's  new  fighter,  the  F-18. 
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Table  1 Reductions  in  Stability  and  Control  Flight 
Testing  Using  Parameter  Identification 


Table  2 Parameter  Identification 
TA-4J  Stability  and  Control  Flight  Tests 


Trim  Airspeed 

<KEAS) 

Number  of 
Maneuvers 

Flight  Time 

165 

2 

1 3 minutes 

200 

2 

250 

2 

300 

2 

350 

2 

400 

2 

V 

MAX 

2 

Total  14 

NOTE:  (1)  Conventional  tests  would  have  required  56- 

maneuvers  and  112  minutes. 


Table  3 Parameter  Identification 
F-14A  Stability  and  Control  Flight  Tests 


Trim  Airspeed 
(KEAS) 

Number  of 
Maneuvers*  * ' 

Flight  Time 

250 

2 

8 minutes 

350 

2 

450 

2 

500 

2 

Total  8 

NOTE:  (1)  Conventional  testing  would  have  required 
32  maneuvers. 
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TabU?  4 Comparison  of  S-3A  Power  Approach 
Parameter  Identification  Results  with  Ongnial 
Contractor/Wind  Tunnel  Data 


Table  6 EA-6B  Phugoid  Parameter  Identification 
Results  in  the  Power  Approach  Configuration(l) 
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Parameter 

Parameter 
Est  imate 

Parameter 
Estimate 
Confidence 
Bound  Hal 

Conventional 

Flight 

Test 

Estimate  * ^ 
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X 
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-0.Z31 
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-0.281 

z 
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0.022 

-1.632 
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00 
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-0.0Z9 

0.002 
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C*a 

*0.271  U> 

- 

- 

c 

X 

q 

0(Z) 

- 

- 

c 

*8 

e 

0UI 

- 

- 

NOTES:  (1) 
(Z) 

(3) 

(4) 


Loading  - 3 ECM  Pods/Z  Tanks 

Wind  tunnel  data  (not  estimated) 
Estimates  are  based  on  flight  test  drag 
polar  data 

Body  fixed  axis  system 


Table  7 


EA-6B  Frequency  and  Damping  Characteristics 
in  the  Power  Approach  Configuration 


Table  5 S-3A  Power  Approach  Modal  Characteristics 


NOTES:  (1)  Loading  - 3 ECM  Pods/2  Tanks 

(2)  Estimate  based  on  combination  of  wind 
tunnel  and  conventional  flight  test  results 


NOTES:  ! Hand  !.chniqu.« 

21  Poaitiv.  tlgn  mdirat**  diT*rg»nt  apiral  mod* 
3>  N.gatir.  »igti  indicates  conrvrgMit  roll  aod» 
141  Clvan  loading 
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Char  art.  rtstK  » 

Conventional 
P.*uit»  1,1 

Wind 

Tuin.l 

Data 

Ai.-frama 

Contractor 
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RrwIU 

spiral  Mod.  Tim* 
r ona'ant  *.< 

Slightly  Dir.rg.nt 
to  Nvutral 

157.73 

17.00 

18.1 

Roil  Mod.  Tim* 
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-0.521 

-0.534 

r inatant 

Dutch  Poll 
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1 .2*7 
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Fraqucnry 

< rad/anr ; 

Dutch  Roll 
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0.102 

0.113 

Damping 
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Wind 

Tunnel 
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Period 

wn 

sp 

( rad/sec \ 

1-071 

1.40b 

1.178 

^sp 

0.630 

0.670 

0.345 

Phugoid 

wn 

P 

( rad/sec ) 

0.148 

0.147(Z) 

0.147121 

cr 

0.092 

0.093(J) 

0.093(Z’ 
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Figure  1 Parameter  Identification  Procedure 
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Figure  2 Iterative  Parameter  Identification  Process 


Figure  3 Root  Locus  Plot  of  the  Roll  Rate  to 
Aileron  Transfer  Function 
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Figure  4 F-14A  Short  Period  Stability  and  Control 

Derivative  Estimates  in  Configuration  Cruise 
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ALTITUDE  — 30,000  FT 

GROSS  WEIGHT  RANGE  - 46,500  - 47,500  LB 
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Figure  5 F-14A  Short  Period  Stability  and  Control 

Characteristics  in  Configuration  Cruise 
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Figure  7 S-3A  Power  Approach  Spoiler  Coefficients 
as  a Function  of  Spoiler  Position  for  Body  Axis 


Figure  6 F-14A  Short  Period  Time  History  Comparison 

in  Configuration  Cruise 
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Figure  9 Comparison  of  S-3A  Power  Approach 
Parameter  Identification  Results  and  Parameter 
Estimates  from  Differential  Engine  Thrust 
Sideslip  Maneuver 
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Figure  1 1 EA-6B  Power  Approach  Short  Period 
Parameter  Identification  Results 
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Figure  10  Comparison  of  S-3A  Power  Approach 
Flight  Test  Results  for  Steady  Heading  Sideslip  and 
Parameter  Identification  Results 
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Figure  12  EA-6B  Power  Approach  Short  Period 
Time  History  Comparison 


Figure  13  EA-6B  Power  Approach  Phugoid 
Time  History  Comparison 
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